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I. INTRODUCTION 

There has been overwhelming evidence that coherent structures such as vortices, stream- 
ers and zonal flows {m = n = 0, where m and n are the poloidal and toroidal modenumbers 
respectively) play a critical role in determining the overall transport in magnetically confined 
plasmas.-"^ Some of these coherent structures, so called streamers, are radially elongated 
structures that cause intermittent, bursty events, which can mediate significant transport of 
heat and particles, for instance, imposing a large heat load on container walls. Zonal flows 
on the other hand may impede transport by shear decorrelation~i2. The Geodesic Acoustic 
Mode (GAM)^"— is the oscillatory counterpart of the zonal flow [m = n = in the potential 
perturbation, m = 1, n = in the perturbations in density, temperature and parallel ve- 
locity) and thus a much weaker effect on turbulence is expected. Nevertheless experimental 
studies suggest that GAMs {n = 0, m = 1) are related to the L-H transition and transport 
barriers. The GAMs are weakly damped by Landau resonances and moreover this damping 
effect is weaker at the edge suggesting that GAMs are more prominent in the region where 
transport barriers are expected.- 

The electron-temperature-gradient (ETG) mode driven by a combination of electron tem- 
perature gradients and fleld line curvature effects is a likely candidate for driving electron 
heat transport.—^— The ETG driven electron heat transport is determined by short scale 
fluctuations that do not influence ion heat transport and is largely unaffected by the large 
scale flows stabilizing ion-temperature-gradient (ITG) modes. 

In this work the flrst demonstration of a high frequency branch of the geodesic acoustic 
mode (GAM) driven by electron temperature gradient (ETG) modes is presented. We 
have utilized a fluid model for the ETG mode based on the Braghinskii equations with non- 
adiabatic ions including impurities and flnite /3 - effects.— »i^ A new saturation mechanism for 
ETG turbulence through the interaction with high frequency GAMs, balanced by Landau 
damping, is found, resulting in a significantly enhanced ETG turbulence saturation level 
compared to the mixing length estimate. 

The remainder of the paper is organized as follows: In Section II the linear ETG mode in- 
cluding the ion impurity dynamics is presented. The linear high frequency GAM is presented 
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and the non-linear effects are discussed in Section III, whereas the saturation mechanism 
for the ETG turbulence is treated in Section IV. The paper is concluded in Section VI. 

II. THE LINEAR ELECTRON TEMPERATURE GRADIENT MODE 

In this section we will describe the preliminaries of the electron-temperature-gradient 
(ETG) mode which is present under the following restrictions on real frequency and wave 
length: fij < a; ~ o;^ << Qe, k^Ci > u > k\\Cf.- Here VLj are the respective cyclotron 
frequencies, pj the Larmor radii and Cj = ^jTjjvfij the thermal velocities. The diamagnetic 
frequency is w^ ~ kopeCe/Ln, k± and k\\ are the perpendicular and the parallel wavevectors. 
The ETG model consists of a combination of an ion and electron fluid dynamics coupled 
through the quasineutrality including finite /3-effects.— "^^ 

A. Ion and impurity dynamics 

In this section, we will start by describing the ion fluid dynamics in the ETG mode 
description, in the limit u > k\\Ce the ions are stationary along the mean magnetic field B 
(where B = Bqc^^) whereas in the limit k_j_Ci » u, k±pi » 1 the ions are unmagnetized. 
We note that the adiabatic ion response follows from the perpendicular ion momentum 
equation by balancing the linear parts of, 

- eniV(j) = TiVrii, (1) 

and we find 

Ui = -r0. (2) 

In this paper we will use the non-adabatic responses in the limits co < k±C[ < k±Ci, r/ = i /^ 
and Qi < u < Qe are fulfilled for the ions and impurities. In the ETG mode description we 
can utilize the ion and impurity continuity and momentum equations of the form. 



"J 



TTljUj- 



dt 
dt 



+ UjV ■ Vj = 0, and (3) 



+ enjV(j) + TjVuj = 0, (4) 
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where j = i for ions and j = I for impurities. Now, we derive the non-adiabatic ion response 
with Tj = Te/Ti and impurity response with with tj = T^/Tf, respectively. We have thus, 



n 



ZT. 



3 



1 

\ I { I )- I \ ¥,. . f. 



(5) 



l-u;7(fcic2) 

Here Tj and n^ are the mean temperature and density of species (j = e,i,I), where fij = 
6n/ni, hj = 6ni/ni and = e(j)/Tf. are the normalized ion density, impurity density and 
potential fluctuations. We also define Zeff = J2 ^k^k/ J^^k^k ~ (^i + z'^hj)/ne. Next we 
present the electron dynamics and the linear dispersion relation. 

B. The electron model 

The electron dynamics for the toroidal ETG mode are governed by the continuity, par- 
allel momentum and energy equations adapted from the Braghinskii's fluid equations. The 
electron equations are analogous to the ion fluid equations used for the toroidal ITG mode, 

OTl 

—p- + V ■ {rieVE + UeV^e) + V ■ {UeVpe + UeV^e) + V ■ (^ef ||e) = (6) 

2^''~dt' "^ '^^'^^^ -Ve + V ■qe = 0. (7) 

Here we used the definitions g^ = —{5pe/2mefie)e\\ x VTg as the diamagnetic heat flux, 
ve is the E X B drift, v^e is the electron diamagnetic drift velocity, vpe is the polarization 
drift velocity, Vtt is the stress tensor drift velocity, and the derivative is defined as d/dt = 
d/dt + peCei X V0 ■ V. A relation between the parallel current density and the parallel 
component of the vector potential (Jy) can be found using Ampere's law, 

Vli|, = -^J||, (8) 

Simplifying and linearizing Eqs. (El [7]and|8]) we find the evolution equations for the electron 
density and the electron temperature in normalized form, 

CjTI (J ~ ~ ~ 

~~d~ ^^at'^ -(! + (! + Ve)Vl) Ve<P - V||VlA|| + 

^n(co«^^^ + si^^|:)(<^-^e-T,) = 0, (9) 

(/3e/2 - Vi)^ + (1 + r7e)(/3e/2)Ve') i|| + V||(0 - ne - Te) = 0, (10) 

(9-5/19 d\ld~, 219-29 

-T. + -e„ cos^-- + sin^- --T. + (,. - -)--0 - --n, = 0. (11) 
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The variables are normalized according to 



),ra,Te) = {Ln/pe){eS(f)/Teo,SnJnQ,STjTeo) 



A\\ = {2CeLn/ (3eCPe)eA\\/TeO, 

Pe = STinTjBl 



(12) 
(13) 
(14) 



Using the Poisson equation in combination with (J5]) we then find 

Tihi/Ue , {Z'^nj/ne)Tj 



He. = - 



uj"^ /k\_c, 



2 ^2 

i 



+ 



a;V(A:ic?) 



+ klXl^ 



(15) 



First we will consider the linear dynamical equations (jHl [lO] and [TTl) and utilizing Eq. (ITj) in 



the same manner as in Refs 



16 and 



171 and we find a semi-local dispersion relation as follows, 



UJ 



/3e 



Ae + ^(1 + Ae) + (1 - e„(l + Ae)) U^ + 



^eh\k±Pe 



k\pl (w - (1 + ?7e)w^)] ( UJ - -e„w^ ) + 

[I + Ae) (w - |e„W^) - (?7e - I) W^ - fwAe 



W 



(f + fciP^)-f(l + ^e 



IW^ 



(16) 



In the following we will use the notation Ag = rj(nj/ne)/(l — w^/Zc^cf) + r/(Ze//n//ne)/(l — 
cj^/fc^Cj) + A^iAfjg. Note that in the limit Tj = Tg, w < k^Ci, k^X^^ < k±pe < 1 and in 
the absence of impurity ions, Ag ~ 1 and the ions follow the Boltzmann relation in the 
standard ETG mode dynamics. Here A^e = v5V(4vm^e^ is the Debye length, the Debye 
shielding effect is important for Xoe/Pe > 1^^. The dispersion relation Eq. flT6|) is analogous 
to the toroidal ion-temperature-gradient mode dispersion relation ecept that the ion quan- 
tities are exchanged to their electron counterparts. In Eq. f lT6|l . the geometrical quantities 
will be determined using a semi-local analysis by assuming an approximate eigenfunction 
while averageing the geometry dependent quantities along the field line. The form of the 
eigenfunction is assumed to be^^. 



^{9) 



^Stt 



;i + cos6') with \9\ < n. 



(17) 
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In the dispersion relation we will replace fcy = (fcy), k± = {k±) and uo = {ujd) by the 
averages defined through the integrals, 



(^'^ > = im Z/"**^^* = *»' (^ + 1^"' - ''■^^ - T'" + ^"') • t^**' 



in /2 5 5 \ 

/»7r 

iV(^) = / rf^^2 ^22) 



Here a = [iq^R (1 + r/g + (1 + //«)) ji^L^) and /3 = 87rno(Te + T^jB'^ is the plasma /3, g is 
the safety factor and s = rq'/q is the magnetic shear. The a-dependent term above (in Eq. 
T8|) represents the effects of Shafranov shift. 



III. MODELING HIGH FREQUENCY GEODESIC ACOUSTIC MODES 

The Geodesic Acoustic Modes are the m = n = 0, kr ^ perturbation of the potential 
field and the n = 0, m = 1, kr ^ perturbation in the density, temperatures and the 
magnetic field perturbations. The high frequency GAM (g, Q) induced by ETG modes {k, u) 
is prevailing under the conditions when the ETG mode real frequence satisfies {Qg > u) > Qi) 
at the scale {k±pi < 1) and the real frequence of the GAM fulfills (fig > fi > fij) at the 
scale {q^ < k^). 



A. Linear Geodesic Acoustic Modes 

We start by deriving the linear electron GAM dispersion relation, by writing the m = 1 
equations for the density, parallel component of the vector potential, temperature and the 
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m = of the electrostatic potential 

- rf-^ + 6„sin^|-0g) - V||Vli(;^ = 0, (23) 

((/3e/2 - Vl)| + (1 + r/e)(/3e/2)V,) i« - V||(n« + f^'J) = 0, (24) 

Q^^eG g^'^eG-U, (25) 

-Vi|0g^-6„sin^|-(n«+T:!;^) = O. (26) 

First we will derive the linear GAM frequency assuming electrostatic GAMs (/3e — ?■ 0) this 
yields a relation between the parallel component of the vector potential and the density and 
electron perturbations using Eq. fl2^ as, 



-Vl^-V„(,lg'+i:f»') = 0. (27) 

The m = 1 component of the electron density can be eliminated by taking a time derivative 
of Eq. fl26|) and using Eq. fl23l) and we get, 



^^^Vi0(°) + ^nv.i^^e^ fe„t;.sin^^ + V,,^ 
ot^ or \ or eno 



PeT^'^l^''^ + ^nv.i^^e- I e„t;.sin^^ + V||^ | ) = 0. (28) 



Here (■ ■ ■ ) is the average over the poloidal angle 9. In the simplest case this leads to the 
dispersion relation. 

Note that the linear electron GAM is purely oscillating analogously to its ion counterpart. 
In previous section we computed the linear dispersion relation for the GAM now we will 
study the non-linear contributions through a modulational instability analysis. 



B. The Non-linear ly Driven Geodesic Acoustic modes 

We will now study the system including the non-linear terms and derive the electron GAM 
growth rate. The non-linear extension to the evolution equations presented previously in 
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Eqs ([9]) - (mD are 

- 1" " ^^1*^ - (^ + (^ + ^^)^^) ^^'^ ~ ^ll^^^ll + 

/ 1 <9 d\ ~ ^ 

en {cose-— + sm6—\ {(j) - Ue - T^) = 

-(/3e/2)[i||,VSi||] + [0,V20], (30) 

('(/3,/2 - Vl)^ + (1 + Ve){Pj2)Vy] i|| + V||(0 - ne-) = -(/3e/2)[0 - He, i||] 

+ (/3e/2)[te,i||] + [0,Vli||], (31) 

-T, + -e„ ^cos^-- + sin^-J -_T, + (r,, - -)--0 - --n, = -[0,Te]. (32) 

In order to find the relevant equations for the electron GAM dynamics we consider the m = 1 
component of Eqs (l30l) - ( l32l) and we find, 



-r.^ + e„sin^|-0g)- V||Vl4^= ([0.,V^0fc])« - (/3e/2)([i||,, VSi||fc])«, 

(33) 
(/3e/2 - Vl)^ + (1 + r/e)(/3e/2)V,) i« - V||(n(j; + f «) = -(/3e/2)([0. - h,k, A\\,])^'^ 
+ (/3e/2)([t,,,i||,])« + ([0fc, Vli||,])«, (34) 

|^S-^|-2 = -(i0.,^e.])«, (35) 

where superscript (1) over the fluctuating quantities denotes the m = 1 poloidal mode 
number and (■ ■ ■ ) is the average over the fast time and spatial scale of the ETG turbulence 
and that non-linear terms associated with parallel dynamics are small since ^ << 1. We 
now study the m = potential perturbations, 

- Vi^4°^ - e„sin^^(n« + T«) = {[4>j.,V'4>k])^'^ - (/3e/2)([i||,, V^i||fc])°. (36) 



We will now neglect the effects of Debye shielding and make use of quasi-neutrality in the 
plasma of the form (n^ = h^ = n) and subtracting Eq. Q from Eq. flTUl) we find, 

Pe^Vi0 + 7^ - env. (cos^-- + siu^-j (0 - n^ - t) - V|| J|| = N,. (37) 

Here we have defined the non-linear term on the RHS as A^i = p^Ce-zx V(^-VV5^(^+^^-V^'-. 
For the GAM we find the (n = 0, m = 1) component of Eq. (Q as, 

dfP 2 9n«_ (,) 

^r~3^r-^2 • ^^^^ 

8 
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This can be written T^, = |fie + A^2^, where the m = 1 component is determined by an 
integral of the convective non-linear term as N2 = — J dtpsCgZ x V0*-'''' ■ VTe ■ This leads 
to a relation between the m = 1 component of the density and temperature fluctuations 
modified by a non-linear term. We continue by taking the m = 1 component of Eq. ( TTOj) 
andm = of Eq. ^, 

^_^_,„„.„„,f!i.<, (39) 

at eriQ or 

P^|V10(°) + e^vM^e^^ (^n« + Nl^) = Nf\ (40) 



Similar to the operations performed to find the linear electron GAM frequency we eliminate 
the m = 1 component of the electron density by taking a time derivative of Eq. ( HOl) this 
yields, 



Note that this will be modified by the effects of the parallel current density {J\\) and the 
non-linear terms, however we see by inspection that on average the term N^ does not 
contribute whereas the N^ non-linearity may drive the GAM unstable. 

We will use the wave kinetic equationi'^'^'^""— to describe the background short scale ETG 
turbulence for {fl,q) < {oJ,k), where the action density N^ = Ek/\ur\ ~ eo|0fcP/i^r- Here 
eol^fcl; is the total energy in the ETG mode with mode number k where eo = t + kj_ + -ttt- 
In describing the large scale plasma flow dynamics it is assumed that there is a sufficient 
spectral gap between the small scale ETG turbulent fluctuations and the large scale GAM 
flow. The electrostatic potential is represented as a sum of fluctuating and mean quantities 

0(x,x,r,t) = $(x,r) + 0(x,t) (42) 

where $(X, T) is the mean flow potential. The coordinates (X, T), (x, t) are the spatial and 
time coordinates for the mean flows and small scale fluctuations, respectively. The wave 
kinetic equation can be written as, 

9 ^r ^ .^ 9 f , p ^\dNk{x,t) d /-. ^\dNk{x,t) 






dt dkx V / dx dx V / dkx 

= 7kNk{x, t) - AuNkix, tf. (43) 
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We will solve Equation (H3|) by assuming a small perturbation (SN/.) driven by a slow 
variation for the GAM compared to the mean {Nko) such that N^ = N^q + SN^. The 
relevant non-linear terms can be approximated in the following form-^, 

{[4>k, Vi0fc]) ^ ql V KkJ-^6Nk{q, Q), (44) 

k '0 

([i||fe, Vii||,]) ^qiy2 Kky\o—5N,{q, fi), (45) 

{[4>k,tk]) ^ -tq.VeT'^—SNkiq^n). (46) 

For all GAMs we have qx > qy, with the following relation between SNk and dNko/dk^, 

^ATfc = -iq^ky<fR—— + ' :, (47) 

where we have used 5ujq = k ■ VEq ~ iikyq^ — kxqy)(p'^ in the wave kinetic equation and 

the definition R = t^ — -. — . Using the results from the wave- kinetic treatment we can 

compute the non-linear contributions to be of the form, 



In order to find the non-linear growth rate of the electron GAM we need to find relations 

(1) _ Cn^x sin 9 



between the variables n^ , Tq and 0^ , 



n 



G ~ r, „2/n''^G' 



n- 



qpn 



(50) 



^G -3^G 3^.2.^:^ — ^(^^_^^^^)i/2^ • (51) 

Using Eqs. ( l50l) and ( 15T1) in the Fourier representation of Eq. ( HTIl resulting in 

fig20g) + e„g.sin^(ng) + 4^)) = ^(0, Vl0)(°), (52) 

and we finally find 



6 " 3 ""^'^ |u;| eo r(r/ - r^^.e)!/^ 

-^A 2V^. -2l^r| p^^fc 

10 
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Eq. ( 153|) is the sought dispersion relation for the electron GAM and we solve it pertubatively 
by assuming i7 = Qq + ^i where Qq is the solution to the linear part c.f. Eq. fl29|) . Now we 
find the perturbation Qi = i'jg which will determine the growth rate of the GAM as, 



7g 



'n 2 2 



Ce/ J-'n 



E 



klp%-i\uj^\NQ 



^ ^n „2 2 



6f2o ■' — ' \u}\ 
5 qlplkyp. 



eo il iVe -VtheY^^ 



'12 a 



2 ^xPe 



/ rb2;fc, 



yHe 



\Ur 



eo n 



dNo 



dK 



12Vi;(r/e-r7the)i/2l-'c|- (54) 

In the last expression we have assumed that the GAM frequency (Qq) can be approximated 
by Qq ^ 2ce/R, i.e. the linear GAM is purely oscillatory as found in Eq. (^U^. The non- 
linearly driven electron GAM is unstable with a growth rate depending on the saturation 
level 10^1 of the ETG mode turbulence. 



IV. SATURATION MECHANISM 



In this section we will estimate a new saturation level for the ETG turbulent electrostatic 
potential (0fc) by using the Landau damping in competition with the non-linear growth rate 
of the GAM in a constant background of ETG mode turbulence, accordingto the well known 
predator-prey models used^^, c.f. Eq. (4) in Ref (jsl) and as well as Ref. (7|), 



o 71 r 



dU, 



G 



dt 



IqUo - IlUg- 



(55) 
(56) 



Here we have represented the ETG mode turbulence as Nj. 






m 



andf/G= (^^sin( 



with the following parameters 7 is the ETG mode growth rate, '-jnl is the non-linear damping 
and 7i is the coupling between the ETG mode and the GAM. The Landau damping rate 
\lL = ItI"^ ) is assumed to be balanced by GAM growth rate Eq. (15^ in stationary state 
^ — )■ and -g^ — )■ 0. In steady state find the saturation level for the ETG turbulent 
intensity as (7, = 7l), 

(57) 



ecpkLn 



A8V2 fLn^^ (L^Y ^T^e 
Pe 



Vthe 



Te Pe ^VtT \qRj \PeJ {qxPefikyPe) 

Note that this saturation level is significantly enhanced compared to the mixing length 
estimate, 



Te Pe 



10. 



(58) 
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Here in this estimation we have used L„ = 0.05, q = 3.0, i? = 4, e„ = 0.025, q^Pe = kyPe = 
0.3 and rje — rjeth ~ 1- Note that the result found using a mixing length estimate is ^-^ ~ 1 
significantly smaller. 

V. CONCLUSION 

In this paper we have presented the first derivation of a high frequency branch of the 
Geodesic Acoustic Mode (GAM). The linear dispersion relation of the high frequency GAM 
showed that the new branch is purely oscillatory with a frequency f] ~ ^. To estimate the 
GAM growth rate, a non-linear treatment based on the wave-kinetic approach was applied. 
The resulting non-linear dispersion relation showed that the high frequency GAM is excited 
in the presence of ETG modes with a growth rate depending on the fiuctuation level of 
the ETG mode turbulence. An analytical expression for the resulting GAM growth rate 
was obtained. To estimate the ETG mode fiuctuation level and GAM growth, a predator- 
prey model was used to describe the coupling between the GAMs and small scale ETG 
turbulence. The stationary point of the coupled system implies that the ETG turbulent 
saturation level (pk can be drastically enhanced by a new saturation mechanism, stemming 
from a balance between the Landau damping and the GAM growth rate. This may result 
in highly elevated particle and electron heat transport, relevant for the edge pedestal region 
of H-mode plasmas. 

The present work was based on a fiuid description of ETG mode turbulence, includ- 
ing finite beta electromagnetic effects and retaining non-adiabatic ions. A more accurate 
treatment based on quasi-linear and nonlinear gyrokinetic simulations is left for future work. 
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